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CN 1 Abstract 

q ' The Green function has a complex dependence upon its underlying domain and dif- 

ferential operator. We briefly review Hadamard's formula for the first variation of the 
Green function due to a perturbation of the domain. We then take a different avenue and 
l/") ' approximate the change in the Green function when the Laplacian is perturbed into a 

number of different operators: Helmholtz, Schrodinger and Laplace-Beltrami. 
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1 Introduction 



> 
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The study of boundary value problems presents a significant computational challenge. In theory 
there are numerous formulas, transforms and methods for solving, say, the Dirichlet problem 
on an ellipse for the Laplace-Beltrami operator. In fact, knowledge of a single function — 
the Green function — would reduce the problem to that of calculating an integral. However, 
the computation of a Green function is prohibitively difficult; considering its central role we 
might turn to methods of approximating it instead. One approach aims to replace the Green 
' function of the problem with that of an easier problem which is, in some sense, nearby. There 

are two ways this heuristic can proceed: variation of the underlying domain and variation of 
the relevant operator. 

The purpose of domain variation is clear; the Green function of a region is nearly impossible 
to compute for all but the nicest of domains. Approximating a given domain by a simpler one 
is a natural way to circumvent the difficulty. The resulting alteration to the Green function 
is very complicated, but a formula due to Hadamard gives the first variation, which can be 
sufficient information for some applications. 

For other applications, we might want to approximate the Green function not under pertur- 
bations of the domain, but rather of the operator. This avenue of thought is less studied and 
for good reason; in many applications the differential operator is considered to be well-known. 
If we turn our attention to inverse problems however, we might pose questions wherein an 
unknown operator has a Green function with some given property. The existence of the under- 
lying operator can be studied by understanding what sorts of variations in the Green function 
can result from perturbing a well-understood operator, such as the Laplacian. 

In this article we briefly review a formal argument developing Hadamard's variation. With 
this in hand we can make a few remarks about what sorts of basic facts can follow from such 
a formula. From here we proceed to consider Helmholtz, Schrodinger and Laplace-Beltrami 
operators which are somehow close to the Laplacian. The resulting variation formulas are new — 
or at least not widely known — and have nice forms open to interpretation, both mathematical 
and physical. 
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2 Definitions and Known Results 



There are numerous definitions and conventions surrounding the Green function of a domain. 
For clarity we take a moment to review a few important properties which will be used in the 
following. Let D C C be a bounded domain. A Green function g : D x D — > [—00,00) is a 
map satisfying three properties: 

1. For every w £ D, g(-, w) is harmonic on D\{w} and bounded outside every neighborhood 
of w. 

2. For every w £ D we have g(w,w) = —00 and g(z,w) — (27r) _1 log|z — w\ + 0(1) as 
z — > w. 

3. For any w £ D we have lim 5(2, tu) = for quasi-every ( £ <9£>; that is, for all £ except 
a set of capacity zero. 

The Green function exists and is unique for any domain whose boundary has positive capacity. 
In particular, we will assume that D has real analytic boundary. Occasionally we wish to 
emphasize the dependence on one argument; hence we sometimes write g w (z) in place of 
g(z,w). We have the following facts: 

1. For any £ £ dD which possesses a barrier, g(z,w) = as z — >• (; having assumed dD real 
analytic, this follows for all (. We can thus extend g to allow one of its arguments to be 
on dD. 

2. The Green function is the solution to the boundary value problem 

{Ag„, =S W in D 
g w = on dD 

It follows that the solution to a more general problem 

(Au = f inD 
1 u = on dD 

is given by integration against the Green function: 

u ( z ) = / 9zfdA. 
Jd 

This property can be taken as the definition of the Green function of D. We can define 
the Green function g* of an operator L analogously: 

J Lg* w =S W in D 
\g* w - on dD 

3. The Green function is subharmonic. 

4. For all z,w £ D the Green function is symmetric: g(z, w) = g(w, z). 

5. Along dD we have 

^(C,«))<fa(C) = (iwKC) 

where u> represents the harmonic measure of D and dg/dn represents the £ derivative in 
the outward normal direction. 
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In what follows, we will need the Poisson- Jensen formula for subharmonic functions. Suppose 
u is subharmonic in a neighborhood of D and u ^ — oo. Then for all z e D 

u(z)= f u(C)dw(z,C)+ / g(z,OAu(OdA(0- 

JdD J D 

Even if Au is a distribution, the formula still holds by interpreting the area integral as a 
distributional pairing. Essentially this result follows from the Green identity 



f ( u— v— — j ds — f (uAv — vAu) dA 

JdD \ on dnj J D 



with v = 



3 Domain Variation of the Green Function 

Before proceeding to variations of operators, we examine Hadamard's formula for the variation 
of the Green function due to change of the domain. Suppose that D is slightly enlarged into a 
new domain D* in such a way so that along dD each point is moved along the outward normal 
direction a small distance. The following first-order variation formula gives a relationship 
between the new Green function g* and the old one. Because this is well-known, we only give 
a formal argument; the original derivation can be found in [3] and a more rigorous treatment 
can be found in [3] or in chapter 15 of [T]. 

Theorem 1 (Hadamard's formula). Let p <= C(dD) be a positive function and suppose that 
for e > each point £ € 3D is moved along the outward normal direction a distance ep(£). 
The Green function g* of the new domain D* satisfies 

g*(z,w) - g(z,w) = -e J ^(z,C)^{w,C)p(C) d8 (0 +o(e). 

Proof. Fix w G D and apply the Poisson- Jensen formula to g^ : 

g*(z,w)= [ 0*(C,ti>)dw(z,O+ / g(z,OAg*(t,w)dA(0. 

JdD JD 

Since Ag^ — S w and g* is symmetric in its arguments, this becomes 

g*{z,w)= g*{w,C)du{zX) + g(z,w). (1) 

JdD 

Denote the outward unit normal by n. For £ G dD, £ + ep(()n £ dD*. Therefore 

g*(w,t + ep(On)=0. 
From here we obtain the first-order approximation 

g*(wX) + e^(wX)p(()+o(e)=0. 
Using this approximation ([T]) becomes 

g*(z, w) - g(z, w) = -e j ^ C)§£(«>, CMC) ds(C) + o(e), (2) 

where we have substituted dg / dn{z,C,) ds{C) for du)(z, £). Finally, we remove the dependence 
on g* by iterating the formula; that is, from ([2]) we have g* — g = o(l) and 

g*(z, w) - g(z, w) = -ef ^(z, [g(w, C) + o(l)} (w, Qp(0 *»(0 + <>{e). 
J dD On On 

This gives the result. □ 
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In traditional notation of the calculus of variations, we can write the aforementioned formula 

as 

g*(z,w) - g(z,w) _ f dg f ^dg i 



5q(z, w) := lim 



We will at times use this notation for sake of clarity. 

Example 2. Consider the unit disk DcC. The Green function for a disk of radius R is 



g(z,w) 



1 

2tt 
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R{z — w) 



R 2 — zw 



Suppose that D is enlarged by uniformly increasing the radius by e — that is, the perturbation 
function p is constant along dD. On one hand, a direct computation shows 



r / \ 9g . , 
og{z,w) = q^[Z,W, 



1 1 



R=l 



2ir 1 — zw\ 



The variational formula gives 



6g(z,w) = - — 
1 

~~ ~2^ 



(i-H 2 )(i-]^| 2 ) 
m | C -<|C-*| 2 ds(i) 
(i-H 2 )(i-N 2 ) 



ao (( - w)(<; - z)(( - z)(( - w) iC 
Since \(\ = 1 we can rewrite ( — z — (1 — (z)/(. The above integral can be written as 

i (i-M 2 )(i-M 2 )C 



Sg(z,w) = / 
2tt« J at 



(C-w)(C-z) (l-C^(i-C^) 



■ dC. 



Note that the second factor in the integral is analytic in a neighborhood of D. The Cauchy 
integral theorem gives 



Sg(z,w) = - 



1 

1 



1 



1 



1 



dD z - w \C - z Q-w 
z{l-\w\ 2 ) w(l-|^| 2 ) 



(l-h| 2 )(l-|z| 2 )C 

(i-a)(i-C^) 



d( + o(e) 



z — w 
1 1 - 



1 — zw 

2 



WZ 



2tt |1 - zw\ 2 ' 
so the two answers agree. 

Before proceeding we make a few remarks. 

1. The outward normal derivative of the Green function is positive since, for instance, it is 
the density of the domain's harmonic measure. Therefore the variation is always negative, 
so we conclude that enlarging a domain decreases the Green function at every point. 

2. Let A : C — > R be a positive smooth function and define the operator L = VAV. If we 
alter the definition of the Green function so that Lg = 5, then we can derive another 
variational formula. From Green's identity 



DD 



. , dv du\ 
A I u- v— I ds 

on on , 



{uLv — vLu) dA, 
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we set v = g z to get an analogue of the Poisson- Jensen formula: 

u(z)= [ u(C)A(C)dw(*,C)+ / g(z,t)Au(QdA(£). 

JdD J D 

From here we obtain the first variation of g. 

g*(z, w) - g(z, w) = e J ^-(z, ()^(w, C)A(CMC) ds{Q + <<e). 

3. In Laplacian growth dynamics, a domain containing the origin grows with outward ve- 
locity given by V(Q = dg/dn((, 0). In these circumstances, an infinitesimal time step dt 
causes a domain variation with ep = V dt. As such, the variation formula yields 

2 

ds(Q. 

A generalization of Laplacian growth dynamics is elliptic growth dynamics, wherein we 
use the Green function for L = VAV and take p(() = X(()dg/dn((, 0). In this case the 
previous formula is modified into 

2 

ds(C). 

4. There is an alternative way of defining the Green function, wherein g(z, w) — — In \z — 
w\ + 0(1) as z — > w. Taking this as a definition leads to a few minor changes. Firstly, 
&9w = —2nS w and the Poisson- Jensen formula reads 

u(z) = [ U (0 du>(z, 0-^1 9(z, 0Au(0 dA(0- 

JdD 2lT JD 

This leads to the variation formula 

g*(z, w) - g(z, w) = ±j ^-(z, ()^-(w, CMC) ds(Q + o(e). 

4 Operator Variation of the Green Function 

Now we turn to situations wherein the domain is fixed but rather the underlying operator is 
somehow close to the Laplacian. We will make repeated use of the following integral operator. 

Definition 1. Given a domain D define the operator T to be integration against the Green 
function of D: 

T4>(z)= [ 4>(Z)9(z,OdA(Z). 

J D 

Lemma 3. The operator T is a bounded linear map from L 2 (D) into C(D), the space of 
continuous functions on D. 

Proof. Notice that T is a left inverse to the Laplacian. By elliptic regularity, T maps L 2 into 
H 2 . But in the plane, Sobolev imbedding implies that H 2 (D) C C(D). □ 



dD on 
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dD on 



A(C)|(0,C) 
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4.1 Helmholtz 



Suppose we'd like to approximate the Green function g* of the Helmholtz operator H = A — a, 
where a is a small constant. In terms of the original Green function g we can derive the 
following. 

Theorem 4. Let D be a bounded domain in C and fix w £ D. The Green function g* of the 
Helmholtz operator A — a satisfies 

9w ~ 9w = aTg w + o(a) 

as a — > ; where the convergence of o(a) is uniform in z for each fixed w. Furthermore, a full 
series expansion is given by 



Proof. The function g* w solves the boundary value problem 

\^9*w =5 w + ag* w in D 
\g* w =0 on dD 

By definition of the classical Green function, this implies that 

g* w (z) = T(5 W (0 + ag* w (0) = g w (z) + aTg* w (0- 

If we define the error function / = g^ — g w the previous equation becomes 

/ = aTg* w = aTg w + aTf. (3) 

It remains to show that Tf — > uniformly as a — > 0. For all a smaller than ||T|| _1 , 

/ = a(I - aT)" 1 ^. (4) 

Since g w € L 2 so is (/ — aT)~ l T g w and we can take a — > to find 

||/|| 2 < \a\\\(I - aT)- l Tg w \\ 2 ^ 0. 

That is, / — > in L? as a — > 0. The lemma implies that Tf — > in C(D), hence uniformly. In 
addition we obtain the full series expansion 



= a " T " 



n=0 



if we use the Neumann series for (/ — aT) 1 in (J3J) ■ □ 



4.2 Schrodinger 

In the previous section we deduced a perturbation formula for a Schrodinger operator with 
constant potential function. With slight adjustments to the proof we can handle a general 
potential function; the operator A — u also possesses a perturbation formula when u is small. 

Theorem 5. Let D be a bounded domain in C and fix w € D. Suppose that p is a smooth 
scalar function defined in a neighborhood of D with corresponding multiplication operator P. 
The Green function g* of the Schrodinger operator A — ep satisfies 

9w ~ 9w = tTPg w + o(e) 



G 



as e — > 0, w/iere the convergence of o(e) is uniform in z for each fixed w. Furthermore, a full 
series expansion is given by 

oo 

9t = E e n (TP) n g w . 

Proof. The first steps of the proof are similar to the Hclmholtz case. The function g^ solves 
the boundary value problem 

=S W + epg* w in D 
g* w = on 3D 

By definition of the classical Green function, this implies that 

g* w = T(6 W + epg* w ) = g w + eTPg* w = g w + eTPg w + eTP(g* w - g w ). (5) 

Define the error function / = — g w , so 

/ = eTPg w + eTPf. 

It remains to show that TPf — > uniformly as e — > 0. 

As in our treatment of the Helmholtz operator, T is a bounded linear map L 2 — > H 2 . Since D 
is bounded, p G L°°(D); thus P maps L 2 — > L 2 . Together this means that TP continuously 
maps L 2 -4- _ff 2 and the rest of proof will proceed analogously as before. For all e sufficiently 
small we find that 

f = e(I-eTP)- 1 TPg w . (6) 

Since g w G L 2 we can take e — > 0: 

||/|| 2 < lellKZ-eTPJ^TP^Ha-^O. 

Thus / in T 2 and TP/ in P 2 C C(D). We conclude that T/ -> uniformly. The 
full series expansion again follows from expanding (1 — eTP) -1 in a Neumann series. □ 



4.3 Laplace— Beltrami 

Next we consider the problem of finding the Green function on a domain for the Laplace- 
Beltrami operator L — VAV, where A is a smooth positive function. If A is close to unity that 
the Green function g* can be approximated by the Green function g of the Laplacian. The 
following result makes this notion precise. 

Theorem 6. Fix w e D and suppose that p is a smooth scalar function in a neighborhood of 
D. Given e > we can define \(z) = 1 + ep(z). Then as e — > the Green function g* for 
L = VAV satisfies 

g*(z,w)-g(z,w)=e f p(0V«/(s,0 • Vg(£,w) dA(0 + o(e), (7) 

JD 

where all derivatives are with respect to £. Furthermore, the error term converges uniformly 
for each fixed w. An alternate formula is also true: 

g*(z,w) - g(z,w) = -eg(z,w) ( P(z) ^ H ) + \ J g z g w ApdA + o(e). (8) 
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Proof. We first give a heuristic argument. Once again, if we begin with the correct integration 
by parts formula we can formally manipulate our way to the perturbation formula above. From 
the Green identity 

/ ^ ( u a v 7T~ ) ^ s — / (uLv — vLu) dA 

J 3D \ on on J J D 

we let u — g z and v — g* w to get 



g z Lg* w dA = / g* w Lg z dA. 

D JD 

Note that Lg* w = S w and Lg z = XS Z + VA ■ Wg z to find 

g(z,w) = \(z)g*(z,w) + / g* w VX-Vg z dA, 

JD 

whence 

g* (z, w) - g(z, w) = -ep(z)g* (z,w)-e / 3* Vp ■ Vg 2 

in 

= -ep(z)g(z, w) - e / g w Vp • Vg 2 rfA + ei? e 

where we have defined 

Re = -p(z)(g*(z,w) - g(z,w)) - I (g* w - g w )Vp ■ V 'g z dA. 



(9) 



Here we end the rigor and claim that R c = o(l). Notice that this is the same gap that appeared 
in the derivation of Hadamard's formula above. Using the claim @ becomes 



g*(z,w) - g(z,w) = -ep(z)g(z,w) - e g w Vp ■ Vg z dA + o(e). (10) 

JD 

We can simplify this via integration by parts. Note that g(£,w) = for £ e dD; therefore 
/ Vp-{g w Vg z )dA = - I P V-(g w Vg z )dA 

JD JD 



= P9w^9z dA- pVg z ■ V.g™ dA 

JD JD 

= -p(z)g(z,w) - / pVg z ■ Vg w dA. 

JD 



Inserting this into (fTU|) gives our first perturbation formula ([7]). To deduce (J8J) we take {10 
and reverse the roles of z and w: 



g*{z, w) - g(z, w) = -ep(w)g(z, w) - e g z Vp ■ Vg w dA + o(e). 

JD 

Adding this half of this equation to half of (fTUj) , we recognize a product rule. 

g*{z, w) - g(z, w) = -eg(z, w) ( E^l±&? j _ 1 J V p • V{g z g w ) dA + o(e). 

Since g(z, £)g{£, w) = on 9D, integrating by parts one more time gives the result. 

We now proceed to derive ([8]) in another fashion, sidestepping the need to consider the error 
term R e . Essentially we will use a change of variables to relate the problem to that of the 
Schodinger operator. 
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Define the functions G* w = g* a/X and u = \~ 1/2 A(A 1 / 2 ). Notice that G* w = on 3D and 

VA = V(VA\/A) = 2VAV\/A. (11) 

With this in mind we compute: 

(A - u)G* w = A(. g ;VA) - ug* w V\ 

= g * w aVa + Vaa 5 ; + 2v g * w ■ vVx - 9 * w aVx 

XAg* w + VA • V 3 ; 



= s w /Vx 



= 8 w /yJ\(w). 



From this we conclude that \J \(w)G* w is the Green function of D for the operator A — u. We 
can use our perturbation formula for this operator if we can determine how u depends upon e. 
We begin with a computation; from equation (jlip we hnd that 



|Wa| 2 = 



VA 



JVAJ; 
4A 



2V\ 
A witl 

AA = A\/AVA = 2\/AAVA + 2|V\/A| 2 = 2Au - 



Now we have what we need to relate u and A without roots. 

jVAj 2 
2A ' 

from which it follows that 

_ AA _ |VA| 2 _ eAp/2 _ e 2 |Vp| 2 /4 
U ~ 2A ~ 4A 2 ~~ 1 + ep ~ (1 + ep) 2 ' 

From here we expand in series to find 

n=Q n=0 

- t(-ir f^+^W 

= — +o(e), 

where the error term is uniformly convergent for all small e. Therefore we have 
yj A(w)A(z)sr*(z, w) = g(z, w) + J (^-+o(e)Jg z g w dA + o(e). 
As a map L°° — > i°° the integral operator 



i-> / 

is continuous, so we deduce that 

o(e)g z g w dA = o(e) 
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uniformly. We conclude that 



y/X(w)X(z)g*(z,w) = g(z,w) + - J ^g z g w ApdA + o(e). 
The last step requires a first order approximation of (A(u;)A(z)) -1 / 2 in e. We find that 



y/X(w)X(z] 



= (l + epH)- 1 / 2 (l + ep(z))- 1 / 2 



= 1- 



ep(z) 



+ o(e) 1- 



ep(w) 



o(e) 



Therefore 

g*(z,w) 



1 - e 



1 -e 



+p(w) 



p{z)+p(w) 



o(e). 



o(e) 



g(z,u)) - eg(z,w) 



p(z)+p(w) 



g(z,w) + -J g z g w ApdA + o(e) 
| y g z g w ApdA + o(e), 



as desired. 

Example 7. Again consider D, the unit disk in C. If we take A(£) = 1 + e|£| 2 then 

1 u) 



□ 



« ff (z,«;)=Re^ / |i! 2 



1 



dA(C). 



£ — to 1 — w£ / \ £ — z 1 — z £ , 
We can evaluate the integral as follows. First we assume that z, w ^ and notice that 



i 



i 



£ — w 1 — w£ / £ — w 1 — iu£ ' 
whence the integrand can be rewritten as 

wz w z 



(£-w)(t-z) (l-z£)(t-w) (1-£D)(£-*) (l-flB)(l-fcf) 



(12) 



We will integrate each of these four terms in turn. Following [2] , we recognize the exponential 
transform of D: 

dA(£) = — irwz\ogEj$(w, z) 



— ttwz log(l — wz) — 2ttwz log \z — w\. 



Here we have used the fact that log(a6) = log a + log b (mod iM.); ultimately we only desire the 
real part of each integral. This handles the first term in (TT2"j) . Next consider 



dA(0, 



R (l-z$(t-w) 



where R = D \ D(w,5) and 6 is a small positive number. Since (1 — z£) 1 is antianayltic 
throughout D we can write 



dA(0 



Z 

w 



£ — w 
log(l - zj) dj 
z JdR €-w 2i ' 



dA(0 
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Considering orientation of the inner circle, we have 



OR 



log(l - zQ d£ 
£-w 2i 



log(l - zQ d£ + w 
2i 



log(l ~ z Q d£ 

z JdD(w.S) £ — HI 



2i 



£ — w 

To evaluate the first of these integrals we remark that z/£ € B for £ € <9B so a series expansion 
gives 



w f log(l — z/£) g?£ 7TW ^ 
z Jon £-w 2i z 



z fc d£ 

fc£ fc (£ — w) 27U 



71 W V ^ 
Z ^ 

0. 



Z fe W J 1 



The dominated convergence theorem implies that 



lim — 

<5->o z 



dD(w,S) 



log(l - gg) d£ _ U) 
2z <5->o z 



£ — w 



log(l — zw — <5ze l0 ) 



d6 



nw 



log(l — zw). 



Thus 



Similarly, 



(1 -*£)(£-«>) 



tM(0 = — log(l -zw). 



(l-£w)(£-z) 
The last term in (Tl2l) can be evaluated as such: 



dA(£) = — log(l - zw). 
w 



f 



1 



(1-^)(1-^) 



dA(0 = - 1 ~[d^ 1 - Z ® 



z 
1 

z 

oo 



= E 

OO 

fc=i 

7T 



l-£iu 
log(l - z/Q d£ 
1 - £w 2i 

, fc^-i+ 1 2? 

(zw) fc_1 

log(l — zw). 



Altogether we have 



Sg(z,w) = —Re 



I _ z w 
zw + — H 



1 

zw 



log(l — zw) — 2zwlog \z — w\ 



(13) 



where we have taken complex conjugates of various terms to combine them better. Note 
that this formula only depends upon rotationally-invariant quantities such as zw, as can be 
expected. 

Next we assume that w = and compute 



s g{Z)0) = Re -L i i 



i 



Re 



4vT 2 



£ - z 1 - z£ 



dA(0 



dA(0. 
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The second term is a harmonic function in £; therefore we have 



Jo 1 - zt, 



= 0. 

6=0 



Given a small S > define Dg = B \ D(z, 6) and use Stokes' theorem to write 

= f ILli^ f i£i 2 g 

JdV>€-Z 2i JdD(z,S) i-zH 

= <k-\ T \z + 6e l9 \ 2 d8 
2 Jo 

^4l-|z| 2 ) 

as (5 — > 0. From this we conclude that 

<53(z,0) 



47T 

We remark the same result is obtained by taking w — > in equation (fl3]) . 
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